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Abstract 

We propose an analytic approach to the problem of collisionless magnetic reconnection formu- 
lated as a process of Alfven eigenmodes' generation and dissipation. Alfven eigenmodes are confined 
by the current sheet in the same way that quantum mechanical waves are confined by the tanh 2 
potential. The dynamical timescale of reconnection is the system scale divided by the eigenvalue 
propagation velocity of the n=l mode. The prediction of the n=l mode shows good agreement 
with the in situ measurement of the reconnection-associated Hall fields. 

PACS numbers: 52.35.Vd, 52.30.Ex, 94.30.cp 
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Magnetic reconnection is recognized as a universal process that converts magnetic field 
energy to the kinetic and thermal energy of the plasma in space and laboratory. The 
physical mechanism of energy conversion and the time scale are the fundamental issues 



of reconnection. Energy conversion in the Sweet-Parker model is done by means 

of magnetic field diffusion into the plasma fluid with the characteristic velocity being the 
diffusion velocity. Petschek introduced a mechanism as a complement to diffusion where the 
change in magnetic field propagates as a slow mode shock 3|. The characteristic velocity 
within Petschek model is the shock propagation velocity, which can be much larger than 
the Sweet-Parker's diffusion velocity. In these models the time scale of reconnection is 
estimated as the system scale divided by the characteristic velocity; the reconnection rate is 
measured by the plasma inflow velocity balanced by the characteristic velocity in a steady 
state. 

However, crucial aspects of reconnection in real plasmas such as onset and temporal 
behavior cannot be resolved in steady theories. The general approach to unsteady recon- 
nection has been to Fourier analyze current sheet dynamics in (u, k) space and search for 



instabilities, e.g. the tearing mode {4]. Unstable reconnection modes grow significantly 
on time scale measured by the growth rate. In this letter we present a novel approach to 
time-dependent collisionless reconnection. Collisionless reconnection is described as the 
generation and dissipation of Alfven eigenmodes. Not only can Alfven eigenmodes grow, 
they can also be damped by transferring wave energy into reconnection ion jets. We solve 
linearized two-fluid equations to find the self-consistent evolution of reconnection following 
an initial perturbation in similar spirit of Landau's method . The dynamical timescale 
of reconnection is the system scale divided by the eigenvalue propagation velocity of the 
n=l mode. Both Fourier analysis of instabilities and this theory are limited in the linear 
regime. 

The experimental motivation for this Alfven-eigenmode approach is the recent in situ 
measurement of reconnect ion- asso ciat ed Hal l fields in collisionless space plasmas 



(see Fig{T]) and laboratory plasmas [ill . Il2j |. Hall fields and current were first introduced 
by Sonnerup Il3fl as a steady structure in the diffusion region. Later numerical studies 

n n 

looked at other various perspectives [U, I15| . In our approach, Hall fields and current are 
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FIG. 1: Schematic diagram of Magnetic Reconnection in the current sheet with Hall fields and 
current included. The Hall fields, the Hall current and the magnetic-field-aligned electric field are 
incorporated into the Alfven eigenmode. Predictions of the Alfven eigenmode are compared with 
measurements from the Polar satellite in Fig. [2j Analysis suggests a phase with opposite Hall 
current and Hall B y preceding the phase in figure. 




FIG. 2: Comparison between the observation by Polar (left column) and predictions of the n=l 
Alfven eigenmode in the Harris sheet (right column). The figure of measurements is from Ref. 7). 
Unit length in x is 100km. 



incorporated into the Alfven eigenmode. Predictions of the n=l mode show good agreement 
with in situ measurements (see FigG]). Hall perturbations are smaller than or, at most, 



B y, y, q. 
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121. This 



comparable to the background as indicated by measurements 
fact implies that a linear theory may suffice to explain the essential physics of collisionless 
reconnection. 



The set of collisionless two-fluid equations is 



C V x B = 4vr(Ji + J E 



(1) 
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cV x E = -d t B, (2) 
d t (n s q s ) + V • J s = 0, (3) 

m s n s d t u s + m s n s (u s • V)u s = n s q s (E + x B) - Vn s T s , (4) 

where the subscript s represents the particle species. Ions and electrons have been assumed 
isotropic and isothermal. We also assume a quasi neutral plasma, rii « n e . The coordinate 
system is depicted in Fig. [H with x along the normal of the current sheet, z parallel to the 
background magnetic field, and y aligned with the background current. d y = is assumed. 



16| with background plasma density n s 



Initially the current sheet is set Harris sheet 
nosech 2 (x/a), background magnetic field B z = .B tanh(— x/a) and the background current 
J yo supported by the ion and electron diamagnetic drifts with velocity u so = 2cT s /q s B a. 
Now we solve the linearized two-fluid equations. The time derivatives of the x component of 
01]), in which d t {n s q s ) is eliminated using ([3]) and d t u sy is eliminated using the y component 
of (HD, are 

v 2 

d t E x = 4ir^-[J ix - p 2 (d xx J ix + d xz J iz ) + ur 2 d tt J ix 

+—(d x J ix + d z J iz ) - n^i—^} + u io d x E y , (5) 
Ui B z 

EyC Tfl e C 2 diE x ,n T ^eo T /«\ 

Jex = n e q e -— H — 2- + p e d xz J ez o z J ez , (6) 

B z m i 'YKV A uj e 

where v\ = B 2 z j ^nnimi is the local Alfven speed, uj s = q s B z /m s c is the local gyrofrequency 



of species s, and p s = \jT s /m s uj 2 is the local gyroradius of species s. In ([6]) we have assumed 
p e d x 1 and uj~ 1 d t <C 1. At the current sheet center where electrons are unmagnetized, the 
above assumption is ambiguous and the electron pressure anisotropy may become relevant 
{17I . fisj ] . In the z component of the electron momentum equation, the pressure gradient 



dominates the inertia effect in high f3 plasmas 19|, [20|. Taking time derivative of this 
equation yields 

d t E z = -(d zz J ez + d zx J ex )T e / {n e q 2 e ) + u eo d z E y . (7) 
Eliminating —qiriiUi B x /c — d z riiTi in the z component of ion momentum equation yields 

minid t u iz = qiUiE^l + Ti/T e ) (8) 



Substitute (JTJ), (jSJ), and ([7]) into time derivative of the y component of ([2]). Rescaling x/a — > x 
and Fourier transforming the achieved equation in z and t (dt ~ — iu;, <9 Z ~ zfc z ) gives 

/ \ / to 2 ^ tanh 2 a; \ a 2 B y S_ 



sech 2 x I \k 2 z Vl sech 2 x / p 2 a + p 2 ao k 2 



where Va = B Q / ^/^n mi is the characteristic Alfven speed, Qi = qiB /rriiC, pi 



y^Ti/rrii^l 2 , p ao = ^T e /miVL 2 and S(x,u,k z ) is the Fourier transform of 

4tt (2n qiC d x J ex mid 2 \ 

S = — — — o zx E y - n qiad zzx u iz + d zx — — o ttz Jex • (10) 

c \ B a sech 2 x Ti + T e J 

In (J9]) we have neglected a; 2 / (f2 2 sech 2 x) compared with (u 2 / k 2 z V\ — tanh 2 xj sech 2 x) in coef- 
ficients of By assuming that B y is of low frequency (uj 2 /Q 2 <C 1) and long parallel wavelength 
(k 2 V^/Q 2 <C 1) as indicated by multi-spacecraft measurements 

as 

. We also neglect terms 

on order of 0{m e /mi) in Qj. Setting B y = ipsechx we turn (Q into 

d xx ip + [X - 1 - (R 2 + A) tanh 2 a# = Ssechx/k 2 z , (11) 
where A = , ^ , -R 



( TTTj) is an inhomogeneous Sturm-Liouville equation. A is an unspecified parameter. The 
weight function is sech 2 x. 

The homogeneous form of (fTTl) is a time-independent Schrodinger equation with cor- 
responding total energy E = A — 1 and a potential well l/tanh 2 x, V = R 2 + A. Only 



bound state solutions exist since E < V. The allowed energy levels (see Ref.|2l|], p. 1653) 
E — V — U V + t — (n + l)] 2 yield eigenvalues of A 



X n = {n + l)n + (1 + 2n) Vi2 2 + 1 + 1. (12) 

The eigenfunctions are 



F(-n, 2JV(X n ) + J - n, |fo n + 1|, ^) 

VVifr) = 7 — \i Lj ^-, (13) 

r v ; (e x + e- x ) b " v ; 

where F is the hypergeometric function, b n = \/V(X n ) — E(X n ) = \/ R 2 + 1. 
Eigenfunctions are real and the first two are ipo( x ) = ( e ~ x + e x )~^ R2+1 and 
■01 (a:) = (e _:c + e x ) _%/ 722+1 (e 2:r — l)/(e 2z + 1). The eigenmode is Alfvenic as indicated 
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by its phase velocity A introduced in (fill) . The Harris sheet confines the Alfven eigenmode 
in the same way a tanh 2 (a;) potential well confines a quantum mechanical wave. From the 



perspective of mode conversion theory [6], Kinetic Alfven Wave (KAW) provides a useful in- 
sight to the eigenmode solution. The KAW dispersion relation is oj 2 /{k 2 v 2 A ) = ^ + kl(p 2 +p 2 ) 



13, [20|. We can heuristically achieve a similar Schrodinger equation with tanh potential 
by replacing k x with —id x in the KAW dispersion relation and employing the x dependence 
of v\, p 2 and p 2 from the Harris sheet. 

With the eigenmode solutions we proceed to calculate B y from equation ( TlTi) as 

By = sechx f G(x\x ]u,k z )sechx ^^jY-^-dx , (14) 

J — oo 

where the Green's function G = Tl,n' l l ) ni. x o)' l Pn{.x) /[(A — A n )$ 2 ] is an infinite series of 
eigenfunctions and $ 2 = ip n (x)ip*(x)sech 2 xdx is the normalization constant (see 



Ref. 



2 1| .chapter 7). Let u = uo r + i^^V ~~ * + ,s = u/i. s is the Laplace transform vari- 
able. We replace the Fourier transform in time with a Laplace transform and treat (|14p as 
an initial value problem. For simplicity the initial condition is set as B y \ t =o = 0, d t By\ t= o = 0. 
The Laplace and Fourier inversion of (fT4"j) is 

V A ^l(x o )sechx 



/oo /*oo pt 
/ / 
-oo J —oo J 



xH[(t - t )VAy/\ l R~ x -\z- z \]S(x , t a , z )dt dz dx , (15) 
where H is the unit step function. B y is in the form of superposition of eigenmodes prop- 



agating in ±z. The phase velocity of the nth mode is Vax/A^R -1 . The sources Ui Z ,J ex 
and E y determine the term S and thus B y . Terms Ji X , J ez , E x and E z , grouped with B y as 
parts of the Alfven eigenmode, can be calculated from the sources using equations (pQ), @, 
(jSJ) and ([7]). The other half of the formulation is the response of the sources to the Alfven 
eigenmode. The field-aligned ion jet Ui Z is calculated from (jSJ); J ex is calculated from (jBJ); 
E y , usually called reconnection electric field, is calculated from 

d xx E y - k 2 z a 2 E y = (4na 2 /c 2 )d t J y , (16) 



where (JTBT) and (fTTJl are Fourier transformed in z, 5 io = a/ c 2 mj / (Anqfrio) , and .Se = 
u io ik z J ez — ui e p 2 ik z d x J ez — uiiJ ix + (uiC 2 /4:7TV A )d t E x . ( JTTj) is the sum of the y component 



of (J4|) multiplied by q s /m s . In obtaining (TlTl) we eliminate J ex using (Q, set n« = n e and 
eliminate d t (n e q e ) using ([3]). Setting £^ = (py/tanhx and inserting ffTTj) into ffT6l) yield 



0**^ + [-£; 2 a 2 + V(z)]p = (47ra 2 /c 2 )5 £ /v / tanhx, (18) 

where V(x) = [2cosh(2a;) — l]csch 2 (2a;) — a 2 sech 2 x/5 2 . V(x) is roughly approximated as 
l/(4x 2 ). Equation (jl8p is an inhomogeneous Sturm-Liouville equation. Two independent 
solutions to its homogeneous form are ip% = ^ l /xH < ^\—ik z ax) and (p 2 = ^/xH^\—ik z ax) if 
k z > 0, and (pi = y/xH^' (—ik z ax) and tp 2 = \fxH^\—ik z ax) if k z < 0. Hq and are 
Hankel functions of first and second kind. We select tpi and ip 2 to ensure a real E y , E y (—k z ) = 
E y (k z ). (pi\ x ^-oo = and v^U^oc = satisfy the left and right boundary conditions, 
respectively. The Green's function of ( TTBl is g(x\x Q ; k z ) = — (/?i(x<)^ 2 (a; > )/A((^i, ip 2 ) (see 



Ref. [2l|. chapter 7), where the smaller (larger) of x and x Q . A(<pi,ip 2 ) = T^i/ft 

is the Wronskian of <pi and (p 2 . Using Green's function, we calculate E y as the response to 
the Alfven eigenmode 

E y = Vtanha; f g(x\x ;k z ) — ^ 7ra c -S E dx a . (19) 
J-oo Vtanhx 

Once we obtain E y by Fourier inversion in z, we calculate B Xl B z and J y through ([1]) and 
(J2J). Now we complete a dynamic formulation of collisionless reconnection. The existence 
of parallel E z indicates the breaking of the "frozen-in" condition. 



The physical meaning of the calculation has two reciprocal parts: The sources (E y , J ex , 
and iiiz) generate Alfven eigenmodes propagating outward in ±z; Meanwhile the Alfven 
eigenmodes excite the sources and dissipate. The eigenmodes-sources coupling evolves 
self-consistently following an initial perturbation. Before the system reaches the phase in 
Fig{TJ we expect an eigenmodes-generation phase in which magnetic energy is converted to 
establish Hall fields and Hall currents. We try to produce this phase with a test perturbation 
E y = \E a \ exp[— (x 2 + z 2 )/l 2 } that is associated with suitable <9 t B in extracting magnetic 
energy and changing magnetic topology globally, I can be as large as the system scale. We 
assume that reconnection ion jets are not established (ui z ~ 0) in this phase. We also assume 
that electrons approximately E x B drift in x (J ex ~ n e q e E y c/ B z ) and avoid any evaluation 
around x = 0. Calculation shows S (jTOj) is quadrupole, sign(S) = —sign(xz). Evaluation 
of ([15]) shows that the n=l mode dominates and sign(B y ) = —sign(S). The contribution 
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from n=l mode is B y = Cipi(x)sechxd z exp[— (z/l) 2 ]t 2 , where sign(ipi(x)) = sign(x) and 
C is a negative constant. In this phase the quadrupole Hall B y and Hall current are 
opposite to those in FigJIJ The Hall current J x fa J ix , like the perpendicular current in 
a KAW, is mainly a modified ion polarization current and associated with the increasing 
of inward Hall E x in the region z/l < 1. Electrons move along the magnetic field to 
track drift ions and keep quasi-neutral, producing the consistent parallel Hall current 
J 2 , J z fa J ez . In the parallel direction force balance is approximately true for electrons, 
n e q e E z fa n e q e u eo B x /c + d z n e T e . The contribution to J Z E Z from n e q e u eo B x /c almost cancels 
out in integrating over x. The density gradient in z created by the drift ions is inward 
at inner current sheet and outward on the outer periphery, relating a total J Z E Z > 0. 
J X E X > and J Z E Z > indicate that the eigenmode stored the converted magnetic energy 
in the form of increasing wave energy. The time scale of this phase is the ion polarization 
drift time, equal to the time of establishing the Hall E x . According to (JHJ) ion experiences 
an total outward force qifiiE z (l + Ti/T e ) in the inner current sheet. The system probably 
transits to the eigenmode-dissipation phase when ion acceleration becomes dominating. 
The ion jets' term can dominate S and produce the Hall quadrupole pattern in Fig{TJ In 
this phase J X E X < and J Z E Z < indicate a decrease in the wave energy transfered to 
the accelerated ion jets and the excitation of secondary E y . In the region x/a ^> 1, the 
excited E y goes as A(k z )\j 'tanh xf2 ~ A(k z ) exp(— |A^|a;) according to f|T9l) . The outward 
ion jet picks up energy qiEJ fa qiE x a. We notice q%E x a ~ (l/2)mjV^ in observations [tIIs!, El- 
Equation (|T5l) explicitly indicates several results. In reconnection S is an odd func- 
tion of x and eliminates all even modes. The amplitude factor l^y/X^R) indicates 
that Alfven eigenmodes are created easier in thinner current sheet and that the n=l 
mode dominates. The step function H suggests that reconnection process can extend 
in ±z at the velocity V^y/Xi/R in this eigenmodes-sources coupling formulation. A 
local change can be communicated with the rest of the system over the dynamical 
timescale LRV^ X x . L is the system scale. The dynamical timescale is often related 
to the time taken for a system to respond to a change and move to another equilibrium state. 

Figj2] presents the comparison between the n=l Alfven eigenmode and the measure- 
ment of Hall fields from Polar satellite 1 71. We model the measured current sheet in 
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Ref. 0] as a Harris sheet with parameters n Q = 8cm -3 , B a = 80nT and a = 150km; 
these numbers are from the observation. T« = 5T e is our estimate for a typical current 
sheet. Independent determination of the amplitude and sign of the Hall fields needs past 
information, which is unavailable. Therefore we take the measured amplitude of B y as 
an input. We compare the x dependence of the Hall fields with the form of the n=l 
mode in the present theory. The perturbation assumption (Sn s /n s < l,B y /B < 1 and 
5B Z /B Q <C 1) is roughly satisfied in this case. In a pure n=l mode B y ~ ipi(x)sechx and 
E x ~ (y^R/ cy/Xy) sinh 2 x[l — (pf/a 2 )d xx + (u io /uJi)d x ]B y . E x is estimated with the absence 
of sources. B y and E x show good agreement with observations in the x dependence. A 
minor difference may result from the deviation of density from n G sech 2 x. In addition, given 
the measured amplitude of B y , the calculated amplitude of E x shows good ag reement with 



data. We also suggest searching n=3 mode signals in lab experiments 
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12|. 



In conclusion, this paper addresses the most fundamental issues of reconnection, namely 
the energy conversion mechanism and the time scale. We propose a new mechanism of 
generating and dissipating Alfven eigenmodes for time-dependent collisionless reconnection. 
The dynamical timescale of reconnection, determined as the system scale L divided by 
the eigenvalue propagation velocity Va\/Xi/R, approaches the Alfven transit time L/Va 
as R —>■ 1. This can be much faster than Sweet-Parker and Petschek reconnection models. 
Notice that the physical meaning of the key result is different in each mechanism. Both the 
diffusion velocity and shock propagation velocity are local outward velocities in =br that 
balance the inflow velocity at a certain interface; the eigenvalue propagation velocity in our 
approach is a velocity in ±z at which a local perturbation communicates globally with the 
rest of system. The growth rate measures the time scale for unstable reconnection modes to 
grow significantly; the dynamical timescale in our approach implies an interval over which 
a new equilibrium is achieved. 



The physics of plasma heating and energization of high energy particles in reconnection 
are not resolved in this paper. Isothermal electrons may be an appropriate approximation 
since u/k z v^ h -C 1 |6]. Isothermal ions can be modified to resolve ion heating through 
Landau damping. A kinetic treatment or test particle method is needed to understand the 
formation of high energy particles. 
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